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AN  APPROXIMATE  SOLUTION  FOR  ROCKET  FLIGHT  WITH 
LINEAR-TANGENT  THRUST  ^I^FUlJ^CONTROL 
By  John  S.  I^pKay1 

Lewis  Research  Center,  NASA,  Cleveland,  Ohio 


ABSTRACT 


This  note  presents  an  approximate  solution  to  the  equations  of  motion 
written  with  linear-tangent  thrust  attitude  control.  The  usefulness  of 
this  solution  is  evaluated  "by  comparison  with  accurate  numerical  integra¬ 
tion.  This  comparison  shows  that  the  solution  gives  accuracies  sufficient 
for  preliminary  calculations  for  satellite  "boost  vehid^ps  and  has  a  four- 
to-one  speed  advantage  over  numerical  integration  of  the  same  accuracy 
level . 


INTRODUCTION 


In  order  to  perform  preliminary  design  calculations  for  satellite¬ 
boosting  systems,  a  reasonably  accurate  estimate  of  the  propellant  re¬ 
quirements  is  necessary.  Unfortunately,  obtaining  such  an  estimate 
usually  requires  the  aid  of  trajectory  calculations.  When  trajectory 
data  are  obtained  by  numerical  integration  of  the  nonlinear  equations  of 
motion,  the  computation  time  often  seems  excessive  for  preliminary  cal¬ 
culations.  This  note,  therefore,  presents  a  closed- form,  approximate 
solution  for  the  equations  of  motion  which  is  about  four  times  faster 
than  stepwise  numerical  integration  of  equal  accuracy. 

Since  the  approximate  solution  presented  (details  given  in  the 
appendix)  has  been  developed  for  use  with  satellite-boosting  systems, 
it  therefore  utilizes  an  approximate  form  of  the  thrust-attitude  schedule 
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that  maximizes  the  mass  in  orbit.  This  schedule,  as  has  been  shown  by- 
Fried  (l)  and  other  authors  (e.g.,  (2),  (3)),  is 

tan  =  a  r  bt 

This  relation,  usually  referred  to  as  the  "linear- tangent"  schedule, 
compares  well  with  results  obtained  by  numerical  integration  of  the 
exact  Euler-Lagrange  differential  equations. 

Theoretically,  the  linear-tangent  schedule  should  be  followed  through 
all  phases  of  flight.  However,  the  large  angles  of  attack  that  can  occur 
with  this  schedule  preclude  its  use  at  low  altitudes.  Consequently,  the 
first  stage  of  powered  flight,  where  aerodynamic  forces  are  concentrated, 
must  follow  some  other  thrust-attitude  schedule,  ordinarily  a  zero-lift 
schedule.  It  is  therefore  reasonable  to  neglect  aerodynamic  forces,  for 
mathematical  simplicity,  during  linear-tangent  flight. 

The  vehicle  burnout  conditions  for  orbital  missions  are  usually 
fixed  (e.g.,  a  circular  orbit  at  a  given  altitude).  For  this  case,  it 
would  be  convenient  if  the  closed  form  solution  would  allow  an  inverse- 
type  calculation  for  the  required  trajectory  parameters  (e.g.,  a,  b,  and 
tf)  from  the  known  end  conditions.  This  is  not  possible  with  the  solu¬ 
tion  presented,  and  further  simplifications  do  not  appear  profitable. 

For  example,  a  somewhat  simpler  solution  given  in  (3),  which  has  about 
l/lO  the  accuracy  of  the  present  solution,  also  fails  to  allow  an  inverse 
computation.  Consequently,  the  trajectory  parameters  must  be  solved  for 
iteratively. 

ACCURACY 

Usually,  selection  of  a  burnout  altitude,  velocity,  and  attitude  is 
necessary  to  specify  an  orbital  mission.  Therefore,  at  least  that  many 
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vehicle  and/or  trajectory  parameters  must  be  variable  to  allow  an  iter¬ 
ative  solution.  One  convenient  set  of  variables  is  a,  b,  and  rn  with 
given  first-stage  burnout  conditions. 

Table  I  shows  the  resulting  values  of  i|Tq,  b,  and  r3  for  a  typical 
three-stage  vehicle  and  a  150-naut  I.cal -mile-orbit  mission.  Also  shown 
are  the  results  for  the  same  case  as  computed  by  accurate  numerical 
integration.  The  errors  shown  in  table  I  have  been  found  acceptable  for 
many  preliminary  design  problems. 

As  will  be  seen  in  the  appendix,  certain  assumptions  will  be  required 
which  allow  the  error  to  increase  with  the  burning  time  of  any  stage.  As 
an  example  of  the  errors  that  can  be  caused  by  these  assumptions,  Fig.  1 
shows  the  error  in  top -stage  weight  ratio  for  three-stage  satellite¬ 
boosting  systems  as  a  function  of  the  total  burning  time  of  the  upper  two 
stages.  The  scatter  in  the  points  shown  is  due  to  other  differences  in 
the  vehicles  chosen,  but  the  essential  feature  to  be  noted  is  that  the 
error  increases  rapidly  after  t  =  400  seconds.  This  corresponds  to 
about  200  seconds  per  stage,  which  does  not  seem  an  unreasonable  limita¬ 
tion  for  typical  stages. 
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where  Rq  is  the  radius  at  the  beginning  of  powered  flight  of  the  stage 
under  consideration. 

Using  the  method  of  variation  of  parameters,  equations  (2)  can  be 
solved  in  the  form 

yf 

xf 

where 


(yo  +  A)  cos  ©tf  +  (yo/©  +  B)  Sin  ajt^  | 
(xq  +  D)  cos  a)tf  +  (xq/oj  +  E)  sin  ©t fj 
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Unfortunately  there  is  no  apparent  way  to  integrate  equations  (4) 
without  some  additional  assumption.  Considering  the  magnitude  of  cu, 
it  can  he  seen  that  it  will  never  exceed  =  1.2398X10"3.  Con¬ 

sequently,  if  t  <  200  seconds,  which  is  a  typical  burning  time  for  one 
stage,  cot  will  always  be  less  than  12°.  For  such  small  angles,  it  is 
reasonable  to  assume 

sin  cut  =  cut 
cos  cut  =  1.0 

Using  equations  (5),  equations  (4)  can  be  integrated  to  give 

A  =  Vj(tG  -  tF  sin  \|r  +  H/b) 

B  =  Vj(F  sin  f  -  G)/cu 

D  =  Vs  (tF  cos  t  +  Cr/b) 

J 


(5) 


E  =  Vj(F  cos  V)/a> 


where 


F  =  In 


G  *  In 


COs(l|ff  -  1(f)  cos 


COs(l|f0  -  l|f)  cos 


cos  Voi 


cos  a 

cos  i|r If 


COS  l(f 

H  =  (cos  Vf  -  cos  ijfo)/cos  %  cos  tf 


Using  the  fact  that  the  method  of  variation  of  parameters  assumes 

•  •  •  • 

A  cos  cut  +  B  sin  cut  =  D  cos  cut  +  E  sin  cut  =  0 

the  velocity  components  can  be  written  as 

yf  =  -cu(y0  +  A)  sin  cutf  +  a>(yo/cu  +  B)  cos  ^f 

=  -cu(xq  +  D)  sin  cutf  +  cu(xq/cu  +  E)  cos  cutf 


J 
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NOMENCLATURE 

a  =*  tan  t^q 

b  =  d(tan  i|r)/dt,  l/sec 

F  =  thrust,  lb 

gc  =  32.174  ft/sec2 

I  =  specific  impulse,  sec 

R  =  radius  from  center  of  Earth  =  -/x2  +  y2,  ft 
r  =  weight  ratio  =  Wq/w^ 

t  =  time  measured  from  ignition  of  i^*1  stage,  sec 

Vj  =  gcI,  ft/sec 

W  =  weight  -  Wq  -  W,  lb 

x  =  distance  along  x-axis  of  Fig.  2,  ft 

y  =  distance  along  y-axis  of  Fig.  2,  ft 

p,  =  gravitational  constant  =  1.4077X101®  cu  ft/sec2 

t|t  =  angle  between  thrust  vector  and  x-axis,  deg 

Subscripts: 

e  =  evaluated  at  surface  of  Earth 
f  =  conditions  at  burnout  of  i^h  stage 
n  =  n^h  stage 

0  =*  conditions  at  t  =  0 

Superscripts : 

( ' )  *  d/dt 

(  )  =  conditions  at  t  =  Wq/W 
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FIGURE  LEGENDS 

Fig.  1.  -  Variation  of  third-stage  weight-ratio  error  with  total,  upper- 
stage  burning  time  for  satellite-boosting  vehicles.  Orbit  altitude, 
150  nautical  miles. 

Fig.  2.  -  Earth- centered,  nonrotating  (inertial)  coordinate  system. 
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TABLE  I.  -  ERROR  IN  FLIGHT  PARAMETERS  FOR 
A  THREE -STAGE  SATELLITE-BOOSTING 
VEHICLE.  ORBIT  ALTITUDE, 

150  NAUTICAL  MILES 


1 

Numerical 

integration 

Approximate  solution 

Error 

•f>  Error 

a 

2.3261208 

-0.0459212 

1.97 

*0 

48.714° 

+0.560° 

1.30 

■ 

0.00323681 

+0.00004885 

1.51 

E-910 
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